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A simultaneous analysis is made of the measured rates of ordinary muon capture (OMC) 
and radiative muon capture (RMC) in liquid hydrogen, using theoretical estimates for 
the relevant atomic capture rates that have been obtained in chiral perturbation theory 
with the use of the most recent values of the coupling constants. We reexamine the basic 
formulas for relating the atomic OMC and RMC rates to the liquid-hydrogen OMC and 
RMC rates, respectively. Although the analysis is significantly influenced by ambiguity 
in the molecular state population, we can demonstrate that, while the OMC data can be 
reproduced, the RMC data can be explained only with unrealistic values of the coupling 
constants; the degree of difficulty becomes even more severe when we try to explain the 
OMC and RMC data simultaneously. 
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1. Introduction 

Ordinary and radiative muon capture (OMC and RMC) on a proton 



H~ +p^n + v», (OMC) (1) 
IT + p -» n + v„ + 7, (RMC) (2) 

are fundamental weak-interaction processes in nuclear physics and a primary source of 
information on gp, the induced pseudoscalar coupling constant of the weak nucleon cur- 
rent, see e.g. jl], [J. The most accurate existing measurements of the OMC and RMC 
rates have been carried out using a liquid hydrogen target, which unfortunately makes 
the analysis of the data sensitive to the molecular transition rates in liquid hydrogen. We 
denote by Au q the OMC rate in liquid hydrogen. The experimental value obtained by 
Bardin et al. [[| is 

A^ = 460±20 [s- 1 ] (OMC). (3) 

As for RMC, Jonkmans et al. [§J measured the absolute photon spectrum for E 1 > 60 
MeV and deduced therefrom the partial RMC branching ratio, R 7 , which is the number of 
RMC events (per stopped muon) producing a photon with E 1 > 60 MeV. The measured 
value of Ry is j§ § 

R* xp = (2.10 ± 0.22) x 10~ 8 (RMC). (4) 

Surprisingly, the value of gp deduced in [|, |5| from the RMC data is ~1.5 times larger 
than the PCAC prediction ||. By contrast, the value of gp deduced in from the OMC 
data is in good agreement with the PCAC prediction. 

On the theoretical side, the early estimation of gp was made using PCAC. Heavy- 
baryon chiral perturbation theory (HB%PT), a low-energy effective theory of QCD, allows 
us to go beyond the PCAC approach, but the results of detailed HB^PT calculations 
up to next-to-next-to-leading order (NNLO) essentially agree with those obtained in the 
PCAC approach. Thus the theoretical framework for estimating gp is robust. The key 
quantities in analyzing OMC and RMC are the atomic rates, A s and A t , where A s (A t ) 
is the capture rate for the hyperfme singlet (triplet) state of the [i-p atom. [] The atomic 
rates for OMC and RMC have also been estimated in the framework of HB^PT [||, [10], 
11], |12 , 13|, |14], 15| . The expressions obtained in HB%PT have been found to be essentially 



in agreement with those of the earlier work [|ll], [T7|, [T8|, [L9], |20| . It has also been confirmed 
that the chiral expansion converges rapidly, rendering estimates of the OMC and RMC 
rates obtained in x?T extremely robust. As for the numerical results, however, the earlier 



estimates of the atomic OMC rates, e.g. 0, [17|]) need to be revised because some values 



of the input parameters (g^, g w N, etc.) used in those estimates are now obsolete. In 
Ref. fllQ"| , we provided updated estimates of A® MC and A° MC based on HB^PT (up to 



1 The rates A s and A t are generic symbols for OMC and RMC. When we need to distinguish OMC 
and RMC, we use the symbols Af t MC and Af t MC . 
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Figure 1: Atomic and molecular states relevant to muon capture in liquid hydrogen; X PPfl 
is the transition rate from the atomic singlet state to the ortho p-fi-p molecular state, and 
X op is that from the ortho to para molecular state. 



NNLO). A notable finding in |H| is that the use of the recent larger value of the Gamow- 



Teller coupling constant, gA, gives a value of A° MC that is significantly larger than the 



older value commonly quoted in the literature, see Refs. QTOJ, fig] . 

To make comparison between theory and experiment, one needs to relate the theoret- 
ically calculated atomic OMC and RMC rates to Au q and R J} respectively. For conve- 
nience, we refer to this relation as the A-L (atom-liquid) formula. Bakalov et al. [21] made 
a detailed study of the A-L formula, and they gave an explicit expression for A[ iq (see 
Eq. (56c) in Ref. pi]). 111 our previous work JK| we analyzed A[ iq using the A-L formula 
of Bakalov et al. and found that the best available estimates of the atomic capture rates 
based on HBxPT would lead to a value of Au q that was significantly larger than A^ p . 
We also reported that, by introducing a molecular state mixing parameter, £, considered 
by Weinberg f22fl , it was possible to reproduce A^ p and R^ xp simultaneously. However, 



the A-L formula of Bakalov et al. does not correspond to the experimental condition of 
OMC; to compare with A^ p , the time sequence of the experimental measurement should 
be considered @. 

In this work we reexamine Au q and i? 7 by incorporating into our analysis the ex- 
perimental conditions as well as the updated estimates of the atomic capture rates. In 
particular, we investigate the influence of ambiguity in the transition rate between the 
molecular states (to be discussed below); we also examine the dependence of the results 
on the values of gp and the molecular parameter £ (see below). 

2. Muonic states in liquid hydrogen 

To evaluate Au q and i? 7 from the calculated atomic OMC and RMC rates, we need 
to know the temporal behavior of the various /i-capture components (capture from the 
atomic states and capture from p-fi-p molecular states). Fig. [I] schematically depicts 
various competing atomic and molecular processes occurring in liquid hydrogen. A muon 
stopped in liquid hydrogen quickly forms a muonic atom (/i-p) in the lowest Bohr state. 
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The atomic hyperfine-triplet state (S=l) decays extremely rapidly to the singlet state 
(S=0), with a transition rate Aio — 1.7 x 10 10 s _1 . In the liquid hydrogen target a muonic 
atom and a hydrogen molecule collide with each other and form a p-\x-p molecule with 
the molecule predominantly in its ortho state. We denote by A PPM the transition rate 
from the atomic singlet state to the ortho p-fi-p molecular state. The ortho p-fi-p state 
further decays to the para p-fi-p molecular state. This rate is denoted by A op . Let N s (t), 
N om (t), and N pm (t) represent the numbers of muons at time t in the atomic singlet, 
ortho-molecular, and para-molecular states, respectively. They satisfy coupled kinetic 
equations, see Eq. (54a) in Ref. pif . To integrate these coupled differential equations, we 



need to know the initial conditions. 

For illustration purposes, let us consider a case in which there is one muon in the singlet 
state at t = 0; i.e., N s (0) = 1 and A om (0) = A pm (0) = 0. We then have 

N s (t) = e" A2t , Ar om (t) = -^ (e -A 3 t_ e -A 2i)) 

A 2 — A3 

at f±\ _ X op \ PPPi P -Mt A pAp P/x -A3* 

pm[ ) (A 3 - A 4 )(A 2 - A 4 ) e (A 2 - A 3 )(A 3 - A 4 ) 6 

_l A pA P p M ^_ Aat ^ 

(A2 — As)(A 2 — A 4 ) 

where A 2 = Ao + \ PPll + A° MC + Af MC , A3 = Ao + A op + A^ IC + A^f c , A4 = Ao + A^ IC + 
Ap^f c . Here A is the muon natural decay rate. A^f c and A°^° are the OMC rates 
in the ortho molecular and para molecular states, respectively; similarly for A^f c , and 
Apmf^-B These rates are given by 

Al = 27o (?Af + \A? ) , A F pm = 2 7P (lAf + hf ) , (6) 



where F stands for "OMC" or "RMC" , and 2 7o = 1.009, 2 7p = 1.143 |T[ . 

At this point we discuss the numerical values of X PPfl and A op . The former shows a 
wide scatter in the literature, ranging from A PPM = (1.89 ± 0.20) x 10 6 s _1 to (2.75 ± 



0.25) x 10 6 s _1 [23]. In this work, for the sake of definiteness, we employ the averaged 
value A PPM = 2.5 x 10 6 s _1 (the main point of our argument is not affected by this choice). 
This value is comparable to the muon decay rate Ao = 0.455 x 10 6 s -1 . As regards 
A op , there is a significant difference between the experimental and theoretical values; 



X e x / = (4.1 ± 1.4) x 10 4 s" 1 as compared with A^ = (7.1 ± 1.2) x 10 4 s' 1 



The dominant state for the OMC and RMC measurements is the ortho molecular state 
as is evident from Eq. (0). In both measurements, data taking starts at t = ti ^ 0, and 
it is essential to incorporate this aspect into the A-L formula (see below). Furthermore, 
in the OMC experiment the time dependence of the population of each state plays an 
important role. 



2 Since Af MC , A^ IC A™ 10 are very small, they can be ignored in the calculation of N Si0miPm (t). In 
evaluating the RMC rate itself, however, we need these capture rates; see Eq. (nOI) below. 
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3. Atom-Liquid (A-L) formula for OMC and RMC 

The discussion so far is common for both OMC and RMC, but we now turn to the 
individual discussion of each case. In the OMC experiment (see Fig. 4 in Ref. @), 
beams arrive at the target on the average in a 3 /zs-long burst with repetition rate 3000 
Hz. The data collection typically starts 1 /is after the end of the 3 /xs-long beam burst, 
and the measurement lasts until 306 /is after the end of the beam burst. As mentioned, 
the cascade processes leading to the [i-p ground state and the transition between the 
atomic hyperfme states are extremely fast. One therefore can safely ignore a time lag 
between the muon arrival time and the time at which the [i-p atomic hyperfine-singlet 
state is formed. To proceed with the consideration of OMC, we assume that the average 
time intervals of Ref. || cited above are actual time intervals. Then, provided all the 
muons arrive at the same time, we can choose with no ambiguity that arrival time as the 
origin of time (t = 0) and let t = U, the starting time for data collection, refer to that 
origin. However, the finite duration (if, = 3 [is) of the beam burst causes uncertainty in 
the value of t — t , to be used in Eq. @; U can be anywhere between 1.0 [is and 4.0 [is. 
To account for this muon pulse duration time, U, we assume for simplicity that the beam 
pulse has a rectangular shape. Then, at time t the average number of residual muons are: 

N^t) = i f b dt'N^t-f). (7) 
£5 Jo 

where N^it) = N s (t) + N om (t) + N pm (t). The OMC experiment || counts the number of 
electrons produced by [i~ — > e~z/ e z/ M , and Au q is deduced from the difference between the 
muon decay rate in liquid hydrogen and that in vacuum; the latter is determined from 
the number of positrons produced in [i + — > e + b , e i7 fl . We use the expression of Ref. || (and 
U = 4 [is) 



JOG ffodNe 



dt 

Judtit-U)^. 



where N e (t) is the averaged number of electrons produced at time t and dN ^ = A iV Ai (t). 
Here we have used the fact the duration of the measuring time (306 [is) is long enough to 
be treated as oo. 

On the other hand, for the RMC experiment j|, |5| , the muons essentially arrive one by 
one and the data taking begins at ti = 365 ns. We therefore can neglect the beam burst 
duration time in the RMC case, and we obtain 

jV 7 (oo)-iV 7 fc) 

^ = — nm ' (9) 

Here N 7 (t) is the number of photons obtained by integrating the photon spectrum over 
the interval, 60 < E 1 < 99 MeV, and the production of photons in RMC is determined 
by 



^ = Ar C N s (t) + AZ^NUt) + A™°N pm (t) , (10) 



5 



9a 9wN 




1.267 13.40 


695 11.9 0.891 x 10~ 3 20.1 x 10~ 3 



Table 1: Coupling constants and the atomic capture rates [s x ] used in the present analysis 



where JV 7 (0) = 0. 

4. Numerical results and discussion 

We first give the numerical values of inputs to be used in what follows. Table [3] presents 
the values of the coupling constants and the atomic capture rates. The OMC and RMC 
rates for the hyperfme-singlet and -triplet states have been calculated in HB%PT up to 



NNLO [CO, |I4| and with the use of the most recent values of the coupling constants 



discussed in 10 



We estimate Au q by using the atomic OMC rates in Table |l|. Besides the A-L formula 
in Eq. (M), we consider two others for the sake of comparison; these two A-L formulas are 



that of Bardin et al. M and that of Bakalov et al. [21] . For the ortho-para transition rate 



we employ either A^f or \ h p . The use of A^ p leads to Au q = 460 s -1 with Eq. @ and 
Aug = 459 s _1 with Bardin et al.'s A-L formula. Those values agree with A e ^ q = 460 ± 20 
s -1 . Meanwhile, if we employ A*p, we obtain Au q ~ 421 s _1 with Eq. ([§]), and Au q ~ 419 
s _1 with Bardin et al.'s formula. Thus, Au q is highly sensitive to X op . On the other hand, 
the use of Bakalov et al.'s A-L formula [21 gives too large a value for Au q regardless of 
whether we use \ e Q x p v or A*p ; A;j,j=532 s _1 for A^ p , and An q =518 s _1 for \ h p . As mentioned 
before, the A-L formula of Bakalov et al, which was adopted in our previous work [fLO] , 



corresponds to the choice of t{ = 0, and this choice does not simulate the experimental 
condition. 

An estimate of i? 7 is obtained from Eq. (^) and the atomic RMC rates in Table ^. With 
the use of A^ p , the calculated value of i? 7 is significantly smaller than R^ xp in Eq. (||); 
R^ xp /R!^ w 1.5. If in Eq. (|S]) we use A*p instead of A^ p , then i? 7 is enhanced by about 9% 
but the increase is not large enough to reconcile Rz 1 with R exv . Thus it is not possible to 
reproduce R exv in the existing theoretical framework with the use of the standard set of 
input parameters. In addition, we remark that our results indicate that the sensitivity of 
i? 7 to X op is comparable to that of Au q . 

As mentioned, the atomic capture rates calculated using a phenomenological relativistic 
tree-level mo del pi are consistent with those of HBxPT p], 0, P, |4f (provided the 
former uses the updated value of Qa and the PCAC value of gp). Therefore, the above 
conclusions are not necessarily unique to HBxPT. However, since HB%PT gives A s and 
At with high precision (primarily because the value of gp is strictly restricted by chiral 
symmetry), it allows us to draw much sharper conclusions than the phenomenological 
approach. 

Next, we discuss the sensitivity of A Uq and i? 7 to possible changes in the values of gp 
and the molecular mixing parameter £. In this discussion we use the phenomenological 
model of Fearing ||18|| , a model which admits the variation of gp within a certain range. 
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As discussed by Weinberg [p2fl , the possible mixing of the ortho molecular p-fi-p spin 3/2 
state and spin 1/2 state, parameterized by £, may change the molecular capture rates to 

Al = £AL(l/2) + (1 - e)AL(3/2), (11) 

where F stands for "OMC" or "RMC" ; Af m (l/2) = Af m [see Eq. (|)] and Af m (3/2) 



27oAf. Although the existing theoretical estimate favors £ ~ 1 [[H], 25], we treat it 



here, as we did in Ref. [10], as a parameter to fit the data. In this phenomenological 



model, with the use of A^ p , we can reproduce R^ xp by adopting either gp = lAgp CAC or 
£ = 0.80. However, with the same value of X op , the OMC data requires gp < 1.2gp CAC or 
£ > 0.95. Therefore, it is impossible to simultaneously fit the OMC and the RMC data 
even by adjusting gp and £. If \ op is taken to be smaller than A^ p = 4.1xl0 6 s -1 , then 
it is not impossible to explain A^ p and Rf* p within the phenomenological model with 
a value of gp larger than that of gp CAC and £ < 0.95. However, we cannot attach too 
much significance to this possibility, since HBxPT constrains the value of gp with high 
accuracy, i.e., there is not much room left for adjusting the value of gp. The result of a 
more precise measurement of X op at TRIUMF |26j will shed much light on this issue. 

Our findings are largely in the nature of reconfirming the conclusions stated in one 
way or another in the literature, but a coherent treatment of OMC and RMC in liquid 
hydrogen as described here is hoped to be useful. Our treatment is characterized by the 
use of the best available atomic capture rates obtained in HB^PT, and by an improved 
A-L formula. Although we have presented examples of simulation of the experimental 
conditions, they are only meant to serve illustrative purposes. Definitive analyses can be 
done only by the people who carried out the relevant experiments. Finally, we remark 
that a precise measurement of the OMC rate in hydrogen gas is planned at PSI P3J . This 



experiment would eliminate the ambiguity of the molecular transition rate discussed in 
this paper and directly test the HB^PT prediction |10], [L^ . 
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